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Abstract 

The heavy quark expansion of Quantum Chromodynamics and the strong 
coupUng flux tube picture of nonperturbative glue are employed to develop 
the phenomenology of hybrid meson decays. The decay mechanism explicitly 
couples gluonic degrees of freedom to the pair produced quarks and hence 
does not obey the well known, but model-dependent, selection rule which 
states that hybrids do not decay to pairs of L = mesons. However, the 
nonperturbative nature of gluonic excitations in the flux tube picture leads to 
a new selection rule: light hybrids do not decay to pairs of identical mesons. 
New features of the model are highlighted and partial widths are presented 
for several low lying hybrid states. 
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I. INTRODUCTION 



Although exphcit nonperturbative gluonic degrees of freedom have been anticipated for 
many years, no clear experimental evidence for their existence has emerged. The search 
for nonperturbative glue, in particular as manifested in hybrid mesons, would be greatly 
facilitated by a rudimentary knowledge of the hybrid spectrum and decay characteristics. 
Although it appears that lattice estimates of light quenched hybrid masses are forthcoming 

hadronic decays remain difficult to calculate on the lattice. Thus one is forced to rely 
on model estimates of the couplings of hybrids to ordinary mesons. Historically, there have 
been two approaches to such estimates. The first assumes that hybrids are predominantly 
quark- ant iquark states with an additional constituent gluon and that decays proceed 
via constituent gluon dissociation 0. The second assumes that hybrids are predominantly 
quark-antiquark states moving on an adiabatic surface generated by an excited "fiux tube" 
configuration of glue Decays then proceed by a phenomenological pair production mech- 
anism (the "^Po model") coupled with a fiux tube overlap 0. 

The authors of Refs. |^ noted the existence of a selection rule: TE hybrids do not decay 
to identical S-wave mesons. A similar observation was made by Isgur, Kokoski, and Paton 
(hereafter referred to as IKP), no hybrids decay to identical S-wave mesonsQ. Although 
this assertion has achieved the status of dogma in hadron phenomenology, we stress that 
it is model-dependent. In particular, in the work of IKP the selection rule follows from 
the assumption that the quark pair production operator is completely decorrelated with the 
gluonic modes in the hybrid (as it is in the successful ^Pq model [^]). This implies that 
the hybrid fiux tube degrees of freedom must annihilate the final state mesonic fiux tubes. 
The result is a spatial fiux tube overlap function which defines the region in which the pair 
creation may occur, has a node along the hybrid quark-antiquark axis, and leads to the 
above mentioned selection rule. However, it is entirely feasible for the hybrid fiux tube to 



^Note that the analogue of TE hybrids does not exist in the flux tube model. 
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be correlated with the pair produced quarks. In this case the argument of IKP breaks down 
and the selection rule need not exist. This is precisely what happens in the model we develop 
here. 

Our approach is based on the successful description of the Dirac structure of confinement 
employed in Ref. |^. The idea is to use the heavy quark expansion of the Coulomb gauge 
QCD Hamiltonian to identify relevant operators. The gluonic portion of these are then 
evaluated using a slightly extended version of the flux tube model of Isgur and Paton 0]. In 
the heavy quark limit the quark and antiquark decouple at leading order. We label these x 
and h where /i^ creates a quark and x creates an antiquark. In this approach the interaction 
contains the kinetic energy term: Hint = — J /i^cr-Dx + H.c, where D = iV + gA. It should 
be noted that this interaction only contains terms which produce or annihilate a QQ pair 
(i.e., gluon production or annihilation by a through-going quark line are not present). This 
implies that contributions to the decay which are higher order in the S-matrix expansion are 
also higher order in 1/Mq. Our model is therefore rather simple: calculate the decay of a 
(flux tube model) hybrid into two ordinary mesons to first order in Hint- The essential new 
feature is that the gluon field operator should be expressed in terms of the nonperturbative 
phonon modes of the flux tube model - as opposed to the traditional expansion in plane 
wave gluons. 

II. FLUX TUBE MODEL OF HYBRID DECAYS 

As discussed above, we shall employ the flux tube model of Isgur and Paton |^ to con- 
struct the hybrid states and the effective decay operator. The model is extracted from the 
strong coupling limit of the QCD lattice Hamiltonian. The Hamiltonian is first split into 
blocks of distinct "topologies" (in reference to possible gauge invariant flux tube configura- 
tions) and then adiabatic and small oscillation approximations of the flux tube dynamics 
are made to arrive at an N-body discrete string-like model Hamiltonian for gluonic degrees 
of freedom. This is meant to be operative at intermediate scales a ~ 6^^/^ where the strong 
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coupling is of order unity. The lattice spacing is denoted by a and there are "beads" (or 
links) evenly spaced between the QQ pair. Diagonalizing the Flux Tube Hamiltonian yields 
phonons, a^A? which are labelled by their color (a), mode number (m), and polarization 
(A) which is transverse to the QQ axis. A hybrid may be built of rimx phonons in the m'th 
mode with polarization A = ±. In particular, hybrid states with a single phonon excitation 
are constructed as 

\H) ~ J rfry,^(r)xS,I),t:A(0,^,-0)^"&J(r/2)4(-r/2)a:^tA,|o). (i) 

Spin and flavor indices have been suppressed and color indices are explicit. The factor Xa a' 
in the hybrid wavefunction projects onto states of good parity and charge conjugation. The 
quantum numbers of these states^ are given by P = i]pc{—)^"^^ and C = ripc{—)^"^^"^^ 
where rjpc = X^i-i = ±1 and N = J2m^{^m+ +^m-)- We shall consider low-lying hybrids 
only so that m = 1 in what follows. 

The definition of the hybrid state makes it clear that decays proceed via the vector 
potential portion of the covariant derivative (recall that we work at leading order in 1/Mq). 
Thus we proceed by writing the vector potential in terms of phonon operators such that 
the expressions for the electric and magnetic fields employed in Ref. ^ are recovered in the 
continuum limit. This yields the following effective decay operator 

i7,„, = '-^Y.f dicosiOTt, h\{irQQ)a ■ e,(fQQ) {a^,, - a^^x) xM^Qq), (2) 

V m,A 

where the e(r) are polarization vectors orthogonal to f. The integral is defined along the 
QQ axis only. The integration over the transverse directions yields the factor which 
may be interpreted as the transverse size of the flux tube. Note that the phonon operators 
represent gluonic excitations which are perpendicular to the QQ axis. Although this appears 



^These expressions differ from Isgur and Paton Q because we have adopted the standard def- 
initions for the polarization vectors and the Wigner rotation matrix, following the Jacob- Wick 
conventions. 
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problematical in traditional perturbation theory, it is required here because, in the adiabatic 
limit, the gluonic field configuration must be defined in terms of the quark configuration 
and therefore the field expansion of the vector potential depends on the quark state under 
consideration. 

The decay amplitude for a hybrid H into mesons A and B is then given by: 



{H\H,r,,\AB) = z^-^^- di J drcos(7rO^^^|^e^^H(r)^^(er)^^,((l-Or)- 

;^MrA('/''^'-'/')XA,AeA(i-)-(o-)] (3) 

where (<t) is the matrix element of the Pauli matrices between quark spin wavef unctions. 
This amplitude should be multiplied by the appropriate fiavor overlap and symmetry factor. 

We note the following general properties of the decay amplitude. The operator is nonzero 
only along the hybrid QQ axis - as follows from the structure of the interaction Hamiltonian. 
Thus qq creation occurs on a line joining the original QQ quarks, smeared over the transverse 
size of the fiux tube. This is in contrast to the model of IKP which has a node along the 
QQ axis. Furthermore the spin operator contracts with the fiux tube phonon polarization 
vector, which is absent in the IKP model. Finally, the decay amplitude vanishes when the 
final mesons are identical for any single-phonon hybrid in an odd mode due to the nodal 
structure in the vector potential. Thus one obtains the selection rule: low-lying hybrids do 
not decay to identical mesons. This subsumes the selection rule of IKP so that none of their 
qualitative conclusions are changed. However we also predict, for example, that hybrids do 
not decay to pairs of identical P-wave mesons. 



III. APPLICATIONS 

All strong hadronic decay calculations depend heavily on the hadronic wavefunctions 
employed. This is especially true here when final states with similar spatial wavefunctions 
are considered. We have taken the simple approach of adopting simple harmonic oscillator 
wavefunctions with width parameters chosen to minimize the eigenenergies of a constituent 
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quark model Hamiltonian . Several of the predicted widths were compared to calculations 
with exact wavefunctions - the results differed by less than 10%, leading us to expect 
that the widths presented below are reasonably reliable. Nevertheless dramatically different 
wavefunction width parameters have been used in the literature (cf., the third reference of 
Ref. 0) so that the reader is cautioned that these results are qualitative (pending more 
detailed analysis p). 

In the following, the normalization is fixed to yield a total width for the 7r(1.8) of 220 
MeV, close to the experimental value of 212 ± 37 MeV for the 7r(1.8). This gives ga^ = 2.56 
GeV~^. Notice that this is consistent with the strong coupling relationship which relates the 
string tension to these parameters, b = g'^Cp/^a'^, if one assumes g ^ 1. Table I presents 
the decay widths for isoscalar, isovector, and ss hybrids with J^'" = 1^+, 0""*", and 1 . 
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TABLES 



TABLE I. Dominant Hybrid Decay Modes. 
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+1=0 w(1.9) 1— Trp P 33 
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KK* 
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mixing angle of 34 degrees has been used for the Ki{\.27) and Ki{\A). 
'^The width of the has been incorporated by convoluting with the appropriate Breit-Wigner 
formula. 

'^The flavor structure of the fj and fj' are taken to be {uu + dd) / \/2 and ss respectively. 
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IV. CONCLUSIONS 



We have used an adiabatic approximation to hybrid structure to study the decay of 
gluonic hybrid mesons. Hybrids may be considered as quarks moving on an adiabatic surface 
generated by a flux tube of nonperturbative glue in an excited state. In the heavy quark 
limit of QCD the decay operator is given by the transverse gluon component of the covariant 
derivative. In keeping with the spirit of the flux tube model of Isgur and Paton, we chose to 
evaluate the gluonic portion of the decay amplitude in terms of phonon degrees of freedom. 
The resulting expression differs significantly from that of Ref. ^ - but surprisingly gives 
similar predictions for the total widths of hybrids. The main difference seems to be in 
substantially smaller D-wave amplitudes predicted here (this is reminiscent of the difference 
between the ^Pq model and the ^5*1 model; see the last paper of Ref. 0). Qualitatively, the 
calculation presented here obeys the well-known spin selection rule, spin zero states do not 
decay to pairs of spin zero states, and a new rule: low lying hybrids do not decay to pairs 
of identical mesons. Decays to mesons with similar spatial wavefunctions are also predicted 
to be suppressed. 

Several candidate hybrid states exist. One of these is the 7r(1800) seen in 7r/o(980), 
7r/o(1300), and K{K7i)s |10| decay channels. Significantly, the channels vrp and KK* are 
suppressed. Although a 3S quarkonium state is expected near 1.8 GeV, the decay charac- 



teristics support the identification of this state as a hybrid (see the discussion in Ref. |jTT 



for a similar analysis). In particular the uss is expected to decay strongly to up while this 
mode is practically zero in our approach. The strength seen in 7r/o(1300) is in accord with 
Table I and is again in conflict with a qq identification of the 7r(1800). The interpretation of 
the 7r/o(980) mode is complicated by the large coupling of the /o(980) to the KK channel. 
Indeed, it is likely that the /o(980) is a KK bound state stablized by the /o(1300) [jl2[. 
In view of this, the large 7r/o(980) mode may be due to strong final state interactions in 
KKq K{K7t). a search for isopartner states fj and r)' which are broad (roughly 330 MeV) 
and which do not couple to pp or will be very instructive. 
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Strong evidence for hybrids would be obtained if the 1 seen at 1.9 GeV at BNL [|T3[ is 



confirmed since these quantum numbers are exotic. This state has been seen in nbi and vr/i 
in accordance with our expectations; however, we also expect a reasonably large KKi{lA) 
mode and suggest that this should be investigated (although this might be complicated by 
strong KK final state interactions). Finally, we consider the prospect for vector hybrid 
production at TJNAF. If one assumes p and uj production by vector meson dominance 
followed by t-channel vr exchange then isospin symmetry implies that the ratio of cross 
sections is a{Cj)/a{f)) = 27. Thus it would be expedient to search for the isoscalar vector 
hybrid. This is an interesting state because all of the normally large S+P channels are 
excluded by quantum numbers and it is therefore very narrow, with a predicted width of 
roughly 60 MeV. The predominant decay mode is vrp; we therefore suggest that the search 
for the UJ be conducted in the tt^tt+tt^ final state. 
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